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Abstract 

Motivated by the idea of imposing paralleling computing on solving 
stochastic differential equations (SDEs), we introduce a new Domain 
Decomposition Scheme to solve forward-backward stochastic differ- 
ential equations (FBSDEs) parallely. We reconstruct the Four Step 
Scheme in [23] with some different conditions and then associate it 
with the idea of Domain Decomposition Methods. We also introduce 
a new technique to prove the convergence of Domain Decomposition 
Methods for systems of quasilinear parabolic equations and use it to 
prove the convergence of our scheme for the FBSDEs. 

1 Introduction 

The theory of forward-backward stochastic differential equations (FBS- 
DEs) is a very active field of research since the first work of Pardoux and 
Peng [28] and Antonelli pQ came out in the early 1990s. These equations ap- 
pear in a large number of application fields in stochastic control and financial 
mathematics. We refer to the monograph [11] , [29] for details, further devel- 
opment and applications. Such systems strongly couple a forward stochastic 
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differential equation with a backward one; and they can be written as a kind 
of stochastic two-point boundary value problems 

dX t = b(t, X t , Y t )dt + a(t, X t , Y t )dW t , 
dY t = -b(t, X t , Y t )dt - a(t, X t , Y t )dW t , (1.1) 
X = x,Y T = g(X T ). 

Together with the theoretical studies on the systems (see [I], [7J, [8], [6], 
[23], [25], [2S], [2B]), finding an efficient numerical scheme for FBSDEs has 
also become an important part of the theory. In order to solve a system of 
FBSDEs, we need to use the "decoupling PDE" technique, based on the so- 
called four step scheme (see [23], [25], [26j). In which, the system of FBSDEs 
is associated with a quasilinear parabolic system of the following type 

f ^ + ^=1^^-+ <V9,b{t,x,9) > +b(t,x,6) = 0, in (0,T)xR», 
\ 9(T,x) = g(x), on W 1 , 

(1.2) 

where 8(t, x) is a vector of m components 6 = (9 1 , . . . , 8 m ), m G N. 
From here, there are two directions to solve FBSDEs. The first trend is to 
solve FBSDEs by using the decoupling technique combining with some prob- 
ability methods to avoid treating the PDEs directly (see [3], [5], [8], [27]). 
The second trend is to solve directly the PDEs. The first paper in this di- 
rection is the one of Douglas, Ma, Protter [10], in which the PDE is treated 
by a finite difference method. Later in 2008, Ma, Shen and Zhao proposed a 
new approach based on the Hermite-spectral Method to treat the PDE (see 
[21]), which is then proved to be much more better than the previous one. 
In this paper, we present a new approach, still based on the second trend, to 
the coupled FDSDEs problem, by combining the classical Four Step Scheme 
with Domain Decomposition Methods or Schwarz Methods, with Waveform 
Relaxation. The idea is to impose Parallel Computing on solving SDEs nu- 
merically. We reconstruct the Four Step Scheme with some new conditions 
and then associate it with Schwarz Waveform Relaxation Methods to par- 
allelize the system of quasilinear parabolic equations fl X . 2 p : System ( II. 2p is 
divided into I subproblems, and each problem is solved seperatedly. The 
scheme is then proved to be well-posed and stable. Up to what we know, 
this is the first attempt trying to apply Domain Decomposition Algorithms 
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to stochastic differential equations. 

In the pioneer work [T7], [TS], [IS], P. L. Lions laid the foundations of 
the modern theory of Schwarz Algorithms. With the development of paral- 
lel computers, the interest in Schwarz Methods have grown rapidly, as these 
methods lead to inherently parallel algorithms. However, the problem of con- 
vergence of Schwarz Methods still remains an open problem up to now. In his 
pioneer work [T7], [TS], [IS], P. L. Lions has proved that the classical Schwarz 
Method for Linear Laplace Equation is in fact equivalent to a sequence of 
projections in a Hilbert space. Moreover, he also observed that the Schwarz 
Sequences of linear elliptic equations is related to Minimum Methods over 
product spaces. This observation was used later by L. Badea in [2] to prove 
the convergence of the classical Schwarz Method for a class of linear elliptic 
equations. 

Later, in [H] and [13], M. Gander-A. Stuart and E. Giladi-H. B. Keller 
applied Schwarz Methods to the 1-dimensional linear advection-diffusion 
equation. Refering to the paper [3], they call Schwarz Methods applied 
to parabolic equations by Schwarz Waveform Relaxation Algorithms. The 
techniques of proving the convergence used in these papers were Laplace 
and Fourier Transforms and some explicit calculations. An extension to the 
nonlinear reaction-diffusion equation in 1-dimension was considered in [13]. 
With the hypothese f'{c) < C in [13], proofs of linear convergence on un- 
bounded time domains, and superlinear convergence on finite time intervals 
were then given in case of n subdomains, based on some explicit computa- 
tions on the linearized equations. Another extension to monotone nonlinear 
PDEs in higher dimension was considered by Lui in [20], [2]], [22]. The main 
idea of the papers is based on the well-known Sub-Super Solutions Method in 
the theory of partial differential equations and the initial guesses are usually 
sub or super solutions of the equations. Recently, an extension to Systems 
of Semilinear Reaction-Diffusion Equations was investigated in [S]. This is 
the first paper trying to apply Schwarz Methods to a system of PDEs in 
1-dimension and the proof of convergence is based strongly on the technique 
introduced in [IB] . 

In order to solve FBSDEs by Schwarz Methods, we encounter the system of 
quasilinear parabolic equations (ll.2|) in n-dimension. We then introduce a 
new technique, which allows us to study the convergence of Schwarz Algo- 
rithms for systems of nonlinear equations in n-dimension. 



3 



2 Forward- backward stochastic differential equa- 
tions 

The structure of this section is as follows: In Section 2.1, we will give the 
definition of forward-backward stochastic differential equations, then state 
some results on the existence and uniqueness of the equations; these results 
will be proved in Section 2.2. 

2.1 Existence and uniqueness results 

Let {W t : t > 0} be a d- dimensional Brownian motion defined on the 
probability space P). We define by {&} the a-field generated by W. 

We suppose that {5t} contains all the null sets of 5 and consider the following 
forward-backward SDEs 

X t = x + J*b(s,X s ,Y s )ds + J*a(s,X s ,Y s )dW s , 

Y t = g(X T ) + ft r b(s,X s ,Y s )ds+ tfa(s,X s ,Y s ,Z s )dW s , 

where t belongs to [0, T}; the processes X, Y, Z take values in R n , R m , R mxd , 
respectively and 6, 6, a, <x, g take values in R n ,M m , R nxd , R mxd and R m , 
respectively. 

Since we are only looking for ordinary adapted solutions of the FBSDEs (12. ip 
(i.e. solutions which are {5t}-adapted and square-integrable, and satisfy 
( 12. ip P-almost surely), we can write ( 12. ip in the following form 

dX t = b(t, X t , Y t )dt + a(t, X t , Y t )dW u 

dY t = -b(t, X t , Y t )dt - a(t, X t , Y t , Z t )dW t , (2.2) 

X = x,Y T = g(X T ). 

Now we state the conditions that we impose on ( 12. ip and (12. 2p : 

(Al) The functions b, b, a, a, g are C 1 -functions with bounded partial 

derivaties; and g is bounded in C 2+5 (R m ) for some 5 in (0, 1). 

(A2) The matrix a satisfies 

\<r(t,x,y)\ < C, 

and 

a{t,x,y)a T {t,x,y) > u{\y\)I,V{t,x,y) e [0,T] x R n x R m , 



4 



where v is a positive continuous function, and C is a positive constant. 
(A3) We impose the following assumptions on a: 
There exists a positive continuous function k such that 

sup{ \z\ : a(t,x,y,z) = 0} < K(\y\),V(t,x,y) G [0,T] x R n x R m . (2.3) 

For fixed (t,x,y), the function z >->■ a(t,x,y,z) is bijective, <7 _1 (£, x, y)(C) is 
contiuous with respect to t, x, y and £; and there exists a continuous function 
A from R to R + , such that 

Iff- 1 ^, x, y)(C) - a- x (t, x, y)(0| < A(|y|)|C - CT , VC, C G K mxn , (2-4) 
where a is a constant lying in [1,2). 

(A4) There exists a positive function rj and a positive constant C such that 
for all (t, x, y) in [0, T] x I" x R m , 

|&(t,x,2/)| < 77(H) 

and 

b(t,x,0) < C. 

(A5) We suppose also that for all k G {1, . . . , m} and for all (t, x, y±, . . . , 
2/fc+i, • • • ,Vm) in [0, T] x R» x R- 1 : 

b k (t, x, yi, ... , y k _ x , 0, j/ib+i, . . . , y m ) = 0, 

and b fc (t, x, 2/1, ... , y fc _i, y fc , 7/ fc+1 , . . . , y m ) is decreasing in y fe . 

Assuming that Y t takes the form 6(t,X t ), P-almost surely, for all t in 
[0, T], by the Ito's formula, we can transform the backward SDE in (12. 2p into 
the following system of PDEs 

f + £";=i < V0, &(*, x, 0) > +&(t, x, 6) = 0, in (0, T) x R», 

#(T,x) = g(x), on R n , 

(2.5) 

where we define a T (t, x, 6*) to be the transposed matrix of a(t, x, 9) and (a*j) 
to be |cr(t, x, 9)a T (t, x, 0). The result on the existence and uniqueness of 
solutions of (12. 2 p then follows 
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Theorem 2.1. (Existence and Uniqueness Theorem) Suppose that Condi- 
tions {Al)- (A4) above are satified, then Equation (\2.2\\ admits a unique 
solution (X, Y, Z) defined as follows: 

The process X is the solution of the following forward SDE 

X t = x+ [ b{s,X s ,6{s,X s ))ds+ [ a(s,X s ,9(s,X s ))ds, (2.6) 
Jo Jo 

where 9 is the unique solution of (\2. 5\i . 

The processes Y t , Z t are then 9{t,X t ) and z(t, X t , 9{t, X t ), V6>(£, X t )), 
where z is a smooth mapping from [0,T]xR n x W 71 x R mxn to R mxd satisfying 

£<r(t, x, y) + a(t, x, y, z{t, x, y, £)) = 0, V(*, x, y, f) G [0, T]xl"x R m x ]R mx ". 

(2.7) 

Remark 2.1. Assumptions {Al) - {A4) is similar to Assumptions {Al) - 
{A4) in W^ . However, in §2.$ , a can vary in [1,2) while in the condition 
(2.12) of {A3) a = 1. {A5) will be used later in the proof of convergence 
for the Parallel Four Steps Domain Decomposition Scheme. 



2.2 Proof of existence and uniqueness results 

This subsection is devoted to the proof of the existence and uniqueness 
results stated in Section [2~T1 We first consider the algebraic equation ( 12. 7\i . 
The following result shows that (12. 7p has a solution. 



Proposition 2.1. Under Assumption {A3), Equation i \2. ?j ) has a unique 
continuous solution z, that satisfies the estimate 

\z{t,x,y,0\ < \{\y\)Mo-{t,x,y)\ a +K{\y\), V(t,x,y,£) G [0, T] xR n xR m xR m 

Proof. Since z(t, x,y,£) is equal to a _1 (t, x,y){£o~{t, x,y)) and a -1 , a are 
continuous, then z(t, x, y,£) is continuous. Finally, Conditions ( 12. 3 j) and 
(j23D give 

\z(t,x,y,£)\ < \a-\t,x,y)(£a(t,x,y))-a-\t,x,y)(p)\ + \a~\t, x, y)(0)\ 
< \{\y\)\Ca{t,x,y)r+\K{\y\)\. 

□ 



6 



Now, the following Proposition states a result on the existence and 
uniqueness of solutions to (I2.5p : 

Proposition 2.2. Suppose that (Al) — (AA) hold. Then the system (Iff. 51) 
admits a unique classical solution 8(t, x), such that 9(t, x), ^0(t, x), V9(t, x), 
A6(t, x) are bounded in C((0, T)xl n ). 

Proof. We first recall a useful result. Let a; be a bounded and smooth enough 
domain of W 1 , we consider the system 



f + EIU ^iagfe+ < V0, b(t, x, 0) > +b(t, x, , 

<p(t, x) = g(t, x) , on [0, T] x du, 

<t>(T,x) = g(T,x), on u. 

(2-8) 

Similarly as in Theorem 7.1, Chapter VII of [16] and Lemma 3.2 of [23J, we 
have the following Lemma 

Lemma 2.1. Suppose that all the functions aij, hi, b are smooth, g is bounded 
in C 1+5 ' 2+5 ([0,T] x u) with 5 belongs to (0,1); and for all (t,x,y) G [0,t] x 
i n xR™ we have 



0, in (0,T) x u, 







vi(\y\)I < ( a ij(t, y)) < V2(\y\)I, 

\b{t,x,y)\<fi(\y\), 
d 



^(t,x,y) 



+ 



dy, 



■Oij(t,x,y) 



<MM), 



for some continuous positive functions Vi{-), ^i ), /"(•)/ an d 

\kt,x,v)\<Ci(l + \y\), 

| <b(t,x,y),y> | < C 2 (l+ \y\ 2 ), 

for some positive constants C\, C<i- Then f lff.ffl admits a unique classical 
solution <p(t,x). Moreover, (f)(t,x), -^(j)(t,x), V </>(£, x), A<p(t,x) are bounded 
in C((0, T) x M n ) by a constant which does not depend on u, and there exists 
a positive number 5' in (0, 1) such that <fi belongs to C 1+s ,2+s ((0,T) x uj). 



(2.9) 
(2.10) 
(2.11) 

(2.12) 
(2.13) 
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Now, we will apply this Lemma to our case. First of all, we verify that 
the conditions of Lemma 12.11 hold. We can see that ( 12. 9 p is a consequence 
of (A2), (I2.10p is a consequence of (A4) and (12. 11 P is a consequence of (Al). 
We only need to prove fl2T2|) and fl2~T3l) . 
Conditions {Al) and (AA) infer that 

\b(t,x,y)\ < \b(t,x,y)-b(t,x,0)\ + \b(t,x,0)\<C l (l + \y\), 

and this implies 

| <b(t,x,y),y> | <C[(l + \y\)\y\ <C 2 (l + \y\ 2 ), 

where C[ is a positive constant. Lemma [2. II then implies that there exists a 
solution (j)(t, x) for all u in IR n bounded and smooth enough. 

By a convergence argument similar as in [23J, we deduce that (12.51) ad- 
mits a unique classical solution 9(t,x), such that 9(t,x), 4^9(t,x), V9(t,x), 
A9(t, x) are bounded in C((0, T) x W 1 ) . ' □ 

We consider the forward SDE on X from (12. 2p . with the assumption 
that Y t can be written under the form 9(t, X t ) 

= b{t, X t , 9(t, X t ))dt + a(t, X t , 9(t, X t ))dW t , 

(2.14) 

x. 

From the Lipschitz condition (Al), we can conclude that (I2.14p has a unique 
solution X, which belongs to L 2 (0, T). We then have the following proof of 
Theorem I2.1[ based on the previous two Propositions. 

Proof. The proof is divided into two steps 
Step 1: (X t ,Y t ,Z t ) is a solution of <\2.2ty . 

We note that the existence and uniqueness of X t has already been dis- 
cussed in the previous paragraph, we only need to prove that (Y t ,Z t ) is a 
solution of the backward SDE in the system (12. 2p . It follows from Ito's 
formula that 




dY 



de ^ dd Li, v -\jr\ ^ d 2 9 



at L — ' oxi ' oxiOXj 

i=l *J=1 



dt (2.15) 



+itZ ai ^ x ^ dw f 



dxi 
i=i 1 
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Equations (I2.15p . Proposition 12.11 and the fact that 9 is a solution of (12. 5p 
lead to 



-b(t, X t , Y t )dt - a(t, X t , Y t , Z t )dW t . 



(2.16) 



Hence, (X t , Y t , Z t ) is a solution of (12. 2h . 
Step 2: (X t ,Y t ,Z t ) is unique. 

Let (X*,Y*,Z*) a solution of (Q, and set Y( to be 6(t,X*), Z[ to be 
z(t,X*, 9(t,X;),V9(t,X*)). It suffices to show that Y{ = Y* and Z[ = Z*. 

We first try to get some estimate for the quantity E(\Y t ' — Y t *\ 2 ). Ito's 
formula gives 

d(Y;-Y t -) = 



1=1 



dxi 



d 2 e 

13 dxidxj 

1,3=1 J 



+ b(t,x;,Y t *) 



+ 



89 



At, xi y;) + a(t, xi y;, z; 



i=i 



dW t , 



dt 



By using Ito's formula for the fc-th component of Y^ — Yf, k G {1, . . . , m}, 
we obtain 



d(Y{ k - y; 

89 k 
~8t 



89 



8 2 9 



i=l 



i,j=t 



b k (t,x;,Y t * 



x2(F/ fe - Y* K )dt + 



89 



J2 X?, 17) + **(«, XI Y t *, z; 



i=i 



(2.17) 



x 



dWt 



+ 



89 



g^At, xt, Y t *) + a% xi y;, zd 



i=l 



dt. 



From the previous equation, the fact that 9 is a solution of (12. 5p implies 

d(Y{ k -Y t * k f (2.18) 
n 89 k 

= 2(>f - Y t * k ){J2 ^iV(t,x;,Y t *) - b\t,x;x)\ + 



i=l 



9 



i,j=l 1 3 

+[b k (t,x;,Y;)-b k (t,x;x)}}dt + 



+2(if - Y* k t ' 



d6 



i=l 



dW t 



1=1 



Since z is a solution of ( 12.7ft . we have that 



d0 



£ x t *, 17) + a fe (t, x*, 17, z; 



i=i 



(2.19) 



i= 1 1 

+ [a k (t, X t , Y*, Z* t ) - a k (t, X*, 1?, Z f ')]. 

Since Yf is equal to Y£, then £(|l^, fc - l^ fc | 2 ) = 0. We can infer from (12381) 
and fl2TT9|) that 



E\Yl k -Y* k \ 2 



(2.20) 



2{y k -Y: k )(£—\b\ S) xiY:) ~b\s,x:x)\ + 



i,j=l 1 3 

+ [b k (s,X:,Y;) -b k (s,X:X)}}ds + 

/T n r\nk 

+ [a k (s,X:,Y;,Z:)-a k (s,X:X,Z' s )}} 2 ds. 

Using (Al), we deduce from ( 12 .20ft that there exists a positive constant Mi 
such that 

E{{Y , k _ y+kft + E j\a k {s,X:.X,Zl) - a k {t.XlX, K)?da 
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< M^Vf - Y : k \(\Y:-Y:\ + \z' s -z:\)ds, (2.21) 



which leads to 



e{\y; - y:\ 2 ) + eJ T \a(s,x:,Y s *, z* s ) - &(s,x* s ,y;, z' s )\ 2 d s 
< M 2 E£\Y;-Y;\(\Y;-Y;\ + \z' s -z:\)ds, (2.22) 

where M 2 is a positive constant. Condition (12. 4p implies that 

\Z' t -Z* t \ = ^a^XlY^ZD-a-^XlY^Z'^ 

< X(\Y t \)\a(t, XI Y t *, Z* t ) - a(t, X?, Y t *, Z' t )\<* 

< M 3 \a(t, XI Y t \ ZD - a(t, X*, Y t *, Z' t )\ a , 

where M3 is a positive constant, since 9 is uniformly bounded . 
Therefore 

e{\y; -y;\ 2 ) + m 4 e \z' s -z*\^ds (2.23) 

< M 2 E \Y' S -Y:\{\Y' S -Y:\ + \Z' S -Z:\)ds 

< m 2 e J* \y: - Y;\ 2 d S + m 2 eJ \y: - y:\\z> s - z* s \d S 

< M 2 E \Yg — Y*\ 2 ds + ~^rE / \Y^-Y s *\^ds 

Jt e 2 ~ a Jt 

+e^M 6 E J \Z' S - Z*\^ds. 

where M4, M5, Mq are positive constants and e is a small constant to be cho- 
sen. Notice that the last inequality comes from Young's inequality. Choosing 
e small enough, we get from (12.231) that 

E{\Y;-Y;\ 2 ) < M 7 J\E(\Y:-Y;\ 2 ) + E(\Y:-Y:\ 2 )^]dsl2.24) 

where M7 is positive constant. 

Secondly, we will prove that E(\Y t ' — Y t *\ 2 ) is bounded and then deduce 
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that E{\Y; - Y t *\ 2 ) is equal to 0. Setting G(t) to be E{\Y( - Y t *\ 2 ), we will 
establish a New Gronwall Inequality, to prove that G(t) is bounded. We have 
the following inequality on G 



Now denoting M 7 f?[G(s) + G(s)^ + l]ds by H(t), we deduce that G(t) is 
less than or equal to H(t)\ and as a result H'(t) is bounded by M 7 [H(t) + 
ff(t)5=5 + 1]. This leads to 



Jo p + p 2 ~ a + 1 

Since (K(H(t)))' is bounded by M 7 , then 

K(H(t)) < M 7 T + K(H(0)) = M 8 . 

In addition, we know that K is increasing, which means that H(t) is less 
than or equal to i^ _1 (M 8 ). We then infer that G(t) is bounded by a postive 
constant Mg. Defining M 9 " 1 G'(t) by P (i), we then deduce that P(t) is positive 
and bounded by 1. This leads to 



where Mio, Mu are positive constants; notice that 1 < a < 2. 

By the classical Gronwall's Lemma, P(t) is equal to 0, which implies that Y( 

is coincided with Y t *. From this, we infer that Z' t is coincided with Z\. 

Now, since every solution (X*, Y* , Z t *) of f !2.2p can be written as (X^,9(t, X*), 
z(t, X* ,9(t, Xt))) and since (I2.14p has only one solution, we can conclude 
that (12. 2p has a unique solution. □ 

3 The Parallel Four Step Domain Decompo- 
sition Scheme 

This section is devoted to the construction of the Parallel Four Step Domain 
Decomposition Scheme and its proof of well-posedness and stability. 




H'(t) 



< M 



We define 





12 



3.1 Definition of the Scheme 

We now define a new Parallel Four Step Domain Decomposition Scheme, 

based on the results obtained in Section 2. 

Step 1: Find a smooth mapping z satisfying (12. 7p . 

Step 2: Choose I to be a constant large enough and consider the domain 
Oi = (—l,l) n . On Oi, consider the following problem instead of ( 12. 5 p 

W + £ l r Wa« 5 fX-+ < V0 l Xt,x,9 l ) > +b(t,x,9 l ) = 0, in (0,T) x O h 
6 l (t,x) = g(x), on (0,T) x dO h 
9 l {T,x) = g(x), on O h 

(3.1) 

The same arguments as in the proof of Proposition 12.21 show that (13. ip has a 
unique classical solution 9 l (t,x), where 9 l (t,x), -^9 l (t,x), V9 l (t,x), A9 l (t,x) 
are bounded. Suppose that 9\t,x) = g(x) on (0,T) x (IR n \C^), then these 
arguments also show that 

lim \ \9 — 9\ |(Loc((o.T)xiR™)) m = 0. 

I— ^oo 

Step 3: Solve the equation (13 .ip iteratively in the following manner 

• Divide Oi into / sub domains 

p=i p=i 

where —l = ai<a 2 <b 1 <---<a I < < bj = I. Denote that 
Si = bi — a i+ i for i belongs to {1, — 1} and Lj = fej — for i 
belongs to {1, ... , I}. 

• Choose a bounded initial guess 9 l in C°°(M n ) at step 0. Associate each 
subdomain Q p with a function 6^ such that 9 l p = 9 l on Q p . 
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• Solve the following p-th subproblem at iteration j^q 

' ^ + EI i=1 < V ^ , ^, t ) > ^ W = °> m(0,T)xQ p , 

e l m {.,.,a p ) = ^_ lj9 _ 1 (.,.,a p ), on (0,T) x {-l,l) n ~\ 
< 9 l Ptg {.,.,b p ) = 9 l p+hq _ 1 {.,.,b p ), on(0,T)x(-/,0 B_1 . 

^(f,x) =s(f,x), on (0,T) x (5O,\((0,T) x (-/, X ({a p } U {6 P }))), 
fl pg (T,x) = <?(!», on H,/)"" 1 x (a p ,6 p ). 

(3.2) 

For the extreme subdomain Qi (resp. f2j), we consider the boundary 
condition 6^ x, ai) = g{t,x) on the left (resp. 9 l T q {t,x,bi) = g{t,x) 
on the right) in (13. 2p . 

• Suppose that we stop at the iteration j^q while solving (13.21) . 

The following two theorems insist that Step 2 of the algorithm is well-posed 
and show that the solutions of the subproblems (13.21) converge to the solution 
of the main problem (13. ip when q tends to infinity. 

Theorem 3.1. (Well-posedness Theorem) Suppose that {Al) — {A5) hold, 
then at each iteration j^q in each subdomain j^p, there exists a unique clas- 
sical solution 9 l pq {t,x) for (jff.ffj) . such that 9 pq {t,x), j^9 l p {t, x), V9 l pq {t,x), 
A9 l (t, x) are bounded and the sequence {9 l pq } p€ {i^..jy jqe ^ is uniformly bounded 
(with respect to p and q) in C((0,T) x 0{). 

Theorem 3.2. Under Assumptions {Al) — {Ah), we have the convergence 

lim sup \\9 p g — 9 l \\(L< x >((o.T)x(-i,i) n - 1 x(a p .b p ))) m = 0. (3.3) 
Q-*°° pe{i,...,l} 

Step 4 : We will continue with the values 9 pq , p G {1, . . . ,1} that we have 
got at the end of step 2. 

• Let 9 l be a function defined on [0, T] x 0\ such that 9 l Jt, x) = 9 l p q {t, x) 
on [0, T] x (fi p \(fi p _ 1 Ufi p+1 )) for p E {2, . . . , J-l}, on [0, T] x (fi p \fip_i) 
for p — I, and on [0, T] x (f2 p \f2 p+1 ) for p — 1. We can choose 9 l g {t, x) 
such that it is Lipschitz, differentiable with respect to x and t in R™ 
and R and 

lim \\9 l q — ^||(i / oo([ ,r]xe> ! )) m = 0- (3-4) 
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• Use 9 q , solve the following forward SDE 

X q t =x + [ b q {s,X q )ds + [ a q {s,X q )dW s , (3.5) 
Jo Jo 

where b q is b(t, x, 9 l q (t, x)) and a q (t, x) is a(t, x, 9 l q (t, x)). 

Using the same arguments as the ones used for (12.141) , we can conclude 

that (13.51) has a unique solution in L 2 (0,T). 

Set Yt l =6%Xf) and Zf = z(t, Xf ,9% Xf),V9% Xf)). The follow- 

ing Theorem says that the sequence (X q, \ Y q '\ Z q ' 1 ) converges to (X t , Y t , Z t ) 
as q and I tend to infinity. 

Theorem 3.3. (Convergence Theorem) Suppose that all the assumptions in 
Section 2.1 hold, then as q and I tend to infinity, {X q, \ Y q '\ Z q ' 1 ) converges 
to the solution (X t , Y t , Z t ) of <\2.1\i in the following sense 

lim lim / E(\Xf' 1 - X t \ 2 )dt = 0, 



fT 

lim lim / E(\Y q ' 1 - Y t \ 2 )dt = 0, 

l— >oo q— >oo Jq 
fT 

lim lim / E{\Zf - Z t \ ~)dt = 0. 

!— >oo q— >oo Jq 



3.2 Proof of Theorem [37L 



First of all, we introduce some useful notations which will be needed for 
the proof. We set 

M = max{||^o|| C ([ 0iT ] xR n), | \g\ |c(r™)}, 

and define p p , q (t, x) to be 9 l pq (T — t,x) for p G {1, . . . , 1} and q G N. 
We can reformulate Systems (I3.2p into the following form 

+ E ■ j= i < V/W b(t, x, p m ) > +b(t, x, p m ) = 0, in (0, T) 

Pp,q{; ;Op) = Pp-l,q-l(; - 5 a p)> 011 ( > T ) X [~h 

< P P ,<?(-, fep) = P P +i, g -i(-, bp), on (0, T) x (-1, 1) 11 - 1 , 

Pp , q (t,x)=g(t,x) 1 on (0,T) x (^((-Z,*)"- 1 X ({a p } U {6 P }))), 
P P , g (0, x) = #(0, x), on (-/, Z)™- 1 x (op, fop). 

(3.6) 
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One can see that ( 13. 6 j) are parabolic systems with the initial condition g. 
Now, we will prove the Theorem by induction. 

At step #1, and in the p-th subdomain, using the same argument as 
in Theorem I2.1[ we can prove that (13. 2p admits a unique classical solution 
9 l pl (t,x), where 9 l pl (t,x), j^9 l pl (t,x), V9 l pl (t,x), A9 l pl (t,x) are bounded. 
Consider the following k-th equation of (13.61) . for k in {1, . . . , m} 

-T + £",-=i < V ^-i' 6 (*> x ' > x ' Pp±) = °' in (°' T ) X fi f 

pJ,i(-> •> a v) = pJ-i,o(-> •> 0> o n (0, T) x (-/, l) n -\ 
i(---,y=4,o(-y. on (0,T) x (-l,l) n -\ 
p k pA (t,x)=g k (t,x), on (0,T) x (00,\(H,O n_1 X (K) U {&,})))> 
^ Pp,i(0,x) =g k (0,x), on H,/) n_1 x (a p ,& p ). 

Using {Ah), we deduce from (13. 7p that 

Q k n q2 k 

—if + E a ^a^+ < Vp^,^^,p p ,i) > +c(*^)^,i = o, 

i,j=l 1 1 

where 

£*(t, 3! , ft ,, 1 )-£*(t,g,^ 1 ,...,p*l 1 ,o,p;+ 1 ,...,py 1 ) 



(3.7) 



c(t, x) 



-Jfe if Pp,l + 0> 

otherwise . 



Since c(t, x) is negative, by applying the maximum principle (see [12J) to this 
equation, we can see that the maximum and minimum of p p l can be obtained 
on the boundaries, for all k in {1, . . . , m}. This means that | \p p> i\ \c([o,T\xn p ) 
is bounded by M , and then M is also an upper bound of | \9 l A \ \c([o,T]xn p )- 
Suppose that up to step qo, the unique classical solution 9 p qo {t, x) exists, 
^.aoC*' 35 )' i^,go( t ' x )' V9 l pqo (t,x), A9 l Piqo (t,x) are bounded, and for all p 
in {1,...,/}, ||^i 9o ||c7([o,T]xnp) is bounded by M . We will show that the 
conclusion is still correct for the step qo + 1. The existence and uniqueness of 
Q pm +i can be infered by using the same argument as in step #1 and Theorem 
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12.11 Now, we consider the following equation, for A; in {1, . . . , m}, 



+ 6 fc (t,x,p Pi g 0+ i) = 0, in (0,T) x 



< 



pLo+iM =9 k (t,x), on (0,T) x (80,\((O,T) X (-1, l)"" 1 x ({a p } U {&,}))), 

k Pi, qo +i(°' x ) = # fc (M), on x (a p ,b p ). 

(3.8) 

Again, by a maximum principle argument applied to Equation (|3.8p . we can 
see that the maximum and mininum of p pgo+ i can only be obtained on the 
boundaries, for all k in {1, . . . , m}. However, we know that ||0p igo +il|c([o,T]xfi p ) 
is bounded by Mq for all p in {1, . . . , /}, from the definition of Mq. We then 
deduce that ||0p i(?o+ i||c([o,T]xfi p ) is bounded by M . 

We then conclude that at each step j^q in each subdomain #p, there exists 
a unique classical solution 9 Ptq (t,x) for (13. 2p and the sequence {# Pig } pe T7-geN 
is uniformly bounded with respect to p and q in C((0, T) x C;) by Mo- This 
concludes the proof. 

3.3 Proof of Theorem I3T21 

We divide the proof into 2 steps 
Step 1: An exponential decay estimate 

Setting e Ptg to be 9 l — 9 l , we deduce the system 

f + E- i= i a;, ^)|^+ < Ve p ,„ 6(t, x, ^) > +c(t, x, 9 l m , 6 l ) = 0, 

in (0,T) x Q p , 

e Pi q(.,.,a p ) = e p _i, g _i(., .,a p ), on (0,T) x (-Z,Z) n-1 , 
e Pl9 (., ., ftp) = e p+ i j(? _i(., bp), on (0, T) x (-/, Z)"" 1 , 
e p<q (t,x) = 0, on (0,T) x (9O,\((0,T) x H,/)™" 1 x ({a p } U {6 P }))), 
e p>q (T,x) = 0, on (-/, /) n_1 x (a p ,b p ), 

(3.9) 
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where 



c(t,x,e l P!Q ,e l ) 



d 2 e l 

dxidxj 



+ < W, [b(t, x, e l M ) - b(t, x, e l )\ > +[b(t, x , e l M ) - h(t, x, e% 



Now, defining e p , q (t, x, y, z) to be e P;? (T — t,x,y,z), we change the system 
into 



( d. 



in (0, T) x Q p , 

e P , q (; ., Op) = e p _i >? _i(., ., Op), on (0, T) x (-/, Z)™" 1 , 
€*,?(•> •, b p ) = e p+ i i9 _i(., ., &„), on (0, T) x (-/, J)"" 1 , 
e Pi9 (t,x) = 0, on (0,T) x (90,\(H, l)^ 1 x ({a p } U {&„}))), 
k e p , g (0,x) = 0, on (-U)™" 1 x (a p ,& p ). 



(3.10) 



We define 



$ p , q (t,x) = ^{e k pA f e?cp{P{x n - u) --ft), 



fc=i 



where /?, w, 7 will be fixed below, and consider the following parabolic oper- 
ator 



£($) = — < V*, 6(t, *, 0) >- £ 



a 2 $ 



/?a, n — . (3.11) 



A direct computation gives 

(3.12) 

m 

= ~ ^ n (*' x ' + ^ a ».» _ /^VnXeJJ 2 exp(/3(x„ - w) - 7) 



fc=i 



+2e: 



det. 



-<V6^,6(M,^)>-E^4 + ^ 
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< E K-7 - ^ x ' 9 U ~ P 2a n,n)(e k p>q ) 2 exp(/?(x n - u) - 7) 



fc=l 

~l~2ep j(? c (t, x, $ Pi <j, $ 

where c fc is the fc-th component of the vector c. 
We consider the following term of (I3.12p 

m 

A = £{("7 - a:. ^) + - /3 2 a„, n )(e^) 2 + 2e k p , q c k (t, x, 9 l m , 9 1 )} 

k=l 
m 



fe=l 

n Q2gl,k 



dec ^ dec a 

1,3=1 •> 

+ < V^ fe , [&(*, x, fljj - &(*, x, 9 l )\ > +[&(*, x, 9 l p J - tyt, x, 9 1 )}}} 

m 
k=l 

+ iVkC^IIA^Hc^, + A^) 2 ||W'||c(^) + N 3 (e k Pig ) 2 } 

m 

< ]£{(~7 - Pb n (t, O + /?a n , n - ^ 2 a n , n )(eJ )g ) 2 + iV 4 (eJ, ? ) 2 }, (3.13) 



fe=i 



where iVi, N2, N%, N4 are constants depending only on the coefficients of 
the system and the bound M Q of 9 l and g in C(M n ). Since all solutions 
of the subproblems {9 l pq } are uniformly bounded, A is negative when 7 is 
large enough and (5 is suitable chosen. This implies that £($ P)9 ) is negative. 
According to the maximum principle, the maximum of <& P)9 can only be 
attained on the boundary of the domain. Which means that the maximum 
of 

m 

k=l 

can only be attained on {0} x W 1 ^ 1 x [a p , b p ], on (0, T) x l)^ 1 x 

({a p } U {&„}))) or on ([0,T] x M n_1 x {a p }) U ([0,T] x W 1 ^ 1 x {6 p }). 
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Since $ M (t, x) is equal to on {0} xM"- 1 x [a p , b p ] and on (0, T) x (dOi\((-l, Z) n_1 x 
({a p } U {& p }))), we have the following cases: 
If 1< p < I, 
in 

x)) 2 exp(/3(x n - w) - 7t) (3.14) 

fc=i 

max > (e* Jt, x)) 2 exp(B(a v — u) — "ft) , 

max V(eJ (t, x)) 2 exp(/?(6 p - w) - jt) } . 

If p= 1 

m 

a;)) 2 exp(/3(x n - u) - jt) (3.15) 

m 

x -M «- 1 x{6 1 }f-' 
fe=l 

If p= 7 

m 

^(e 7 fc 9 (t,x)) 2 exp(/?(x n - w) - -yt) (3.16) 

k=l 

m 

< max > (e k r Jt, x)) 2 exp(B(aj — u) — jt). 

~ (t,x)e[0,T]xl-l,l]"^x{aj} ^ 

Step 2: Proof of the Convergence 

Step 2.1: Estimate of the right boudaries of the sub-domains. 

For x in [— l,l] n , we denote x by (X,x n ), where X G [— M] n_1 an d x n £ 
[—/,/]. Moreover, we define 



fc=i 

< 



E Q = max 

P e{i,...,/} I (t,; 



max y^(ep (t,x)) 2 exp(-7t) I . 
:,x)e[o,T] x [-/,/]» ^ v MV ^ v J 

Consider the 7-th domain, at the g-th step, we can see that (13.161) infers 

m 

^{e k Ijq {t, X, x n )) 2 exp(/3(x n - a 7 ) - 7*) 



fc=i 
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< ,„ , yi( e L( t ' x ' a ^) 2ex p(-^)' 

0,2lx -U ^ 



fe=i 



where w is replaced by aj. 

Replacing x n by in the previous inequality, we obtain 

in 

X, bj.,)) 2 exp^&j-! - a z ) - 7 

k=l 

m 

< max V(eJ (t, X, a 7 )) 2 exp(- 7 t). 

fc=i 

Since Jt, is equal to ef_ x X, then 

m 

^(et lj(?+1 (t,X,6 / _ 1 )) 2 exp(/3(^ 1 - a x ) - jt) 
fc=i 

m 

< max V(eJ (t,X,a/)) 2 exp(-7t). 



k=l 



We define /?i to be and let /3 in this case be then if we choose 7 large, 
7 — /3 2 is large, the inequality becomes 



^( e ^ w (t,X,6 / _ 1 )) 2 exp(-7t) 



fc=i 

m 

< exp(- ( 5i5/_i) max VVeJ (t,X, a 7 )) 2 exp(-7£). 

fc=l 

We deduce that 

m 

^(et li9+1 (t,X,6 / _ 1 )) 2 exp(- 7 t) < exp(-faSi-i)E r (3.17) 
fc=i 

Moreover, on the (J — l)-th domain, at the (q + l)-th step, (I3.14p leads 

to 

m 

J^( e i-i,g+i(*>^ a; n)) 2 exp(/3(a; n - o/_i) -7*) 

k=l 
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< max < max > (eL „,i (£, X, a/_i)) 2 exp(— 7t), 

I (t,X)G[0,T]x [-/,/]«-! ^ ^.Si-J-V 
m 

(tfleloflxH,!]"- 1 ^ 
where a; is replaced by a/_i. 

Since ej_ x „ +1 (£, X, 61-2) is equal to ejL 2g+2 (£, X, 61-2), then 

m 

X)( e *-3,ff+3(*' X ' & /-2)) 2 exp(/3(6 / _ 2 - fl/ _i) - 7*) 



fc=i 



< max 



{m 
max y^( e /-i ^ ai-i)) 2 exp(~7£), 
fc=l 



max y^( e /-i 0+1 (*) fo /-i)) 2 exp(/3L/_i - 7*) 

k=l 

thus 

m 

$>*-W^ X, ft/_ 2 )) 2 exp(^5' J _2 - 70 
fe=i 

max y^( e /-i 0+1 (*> ^ a/-i)) 2 exp(-Tt), 
fc=i 

m ~\ 

max y.( e /-i 0+1 (*) x , bi-1)) 2 exp(/3L 7 _i - 7*) > . 

fc=i j 

Combining this with ( 13.1 7p and the fact that 

m 

, ^ ,„ , yi( e /-l,</+l(*' X ' a J-l)) 2eX P( - 7*) < Eq+U 

(t,X)G[0,T]x[-«,i] n ^ 1 ,y 
fc=l 

we obtain that 

m 

y^(ej_ 2 q+2 (t, X, 6/- 2 )) 2 exp(/3S'/_2-7 i ) < max{£q exp^Lj-i-^i^-i), E q+1 }. 
fc=i 
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Thus 



22(e*_ 2 q+2 (t,X, 6/„ 2 )) 2 exp(-7t) < max{E q exp(/3(Lj_i - S/_ 2 ) - P1S1-1), 



k=l 



-E g+ iexp(-/3S/_ 2 )}. 



Denning /3 2 to be and choosing /3 to be /3 2 such that 

+ Sj-2) + ^5/.! = /3 2 S 7 _ 2 , 

we infer 

m 

^(et 2 , g+2 (t,^,&/- 2 )) 2 exp(- 7 t) < max{ J B fc , J B fc+1 }exp(-/3 2> S / _ 2 ). (3.18) 

k=l 

Using the same techniques as the ones that we use to achive (13.171) and 
(13.181) . we can prove that 

m 

^2(e k j_ j q+j (t, X, &/_j)) 2 exp(-7t) < max{E k , . . . , E k+j _x} exp(-/3j £/_,), 

k=l 

(3.19) 

where 

& = A^...y^, J€{2,. ..,/-!}. 
-w-i ^i-j+i 

Step 2.2: Estimate of the left boudaries of the sub-domains 
Consider the 1-th domain, at the k-th step. Then (13. 1 51) infers that 



k=l 

< 



^2{e k 1>q {t, X, x n )f exp(/3(x n - 60 - 7 t) 
in 

m ^ , 1 ^( e L(^ x ' 6 i)) 2ex p(-^)' 

),T]x[-Z,Z] n - 1 z — ' ' y 



(t,X)G[0,T].. 1 , ; | 



we notice here that u is replaced by b\. 
Replace x n by a 2 , we obtain that 



^(e^(t,X,a 2 )) 2 exp(/3(a 2 - b x ) - jt) 



k=l 

23 



0,T x -U ^— ' ,y 



(t,X)G[0,T]xH,/r 

fc=l 



Since e*_(t, X, 02) is equal to e2 j(?+1 (t, X, g^), 
X>W*> X ' a 2 )) 2 exp(- 7 t) 



k=l 

< 



m , a ^ , 1 y]( e i,g( t '^ : ^i)) 2ex P(-7^)exp(/3»S 1 ), 



(t,x)e[o,T]x[-M] 

fe=i 



We define /3[ to be — A/2 an d l e t be /3( in this case. If we choose 7 large, 
7 — P 2 is large. The inequality becomes 



X]( e t+i( t ' X ' a 2)) 2ex P(-70 



fc=i 



< exp(-#£i) max J>* g (t, X, &0) 2 exp(- 7 t). 

fc=i 

We deduce 

^( e ^ +1 (t,X,a 2 )) 2 exp(- 7 t) < exp(-/3^ 1 ) J E; 9 . (3.20) 

k=l 

Moreover, on the 2-th domain, at the (q + l)-th step, (I3.14p leads to 

m 

^(eiq+iiti x i x n)) 2 exp((3(x n - a 2 ) - -ft) 

k=l 

(m 
max >(eo . , 1 (t, X, a 2 )) 2 exp(— 7 £), 



fe=l 



m ? x , 1 y]( e L+i(^ x ' 62 )) 2ex p(/ 3 ( 6 2 - ° 2 ) - t*) > , 



(t,X)e[0,T]x ; | 



notice that u is replaced by a 2 . 

Since e2 j9+1 (t, X, a 3 ) is equal to e 3q+2 (t,X,a 3 ), then 

m 

X^( e 3.g+2(^ X ' a 3 )) 2 exp(/3(a 3 - a 2 ) -7*) 
fc=i 
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< max 



{m 
max y^(e% a+1 (t,X, a 2 )) 2 exp(-'yt) 
(i,x)e[o,T]x[-z,z]»-i ^ ' q+ " v ' 

k=l 



, „ r 1 ^ i 5^( e L+i( t ' X ' fo 2)) 2 exp(^(6 2 -a 2 ) -7t) 
(t,x)e[o,r] x [-y]"- 1 ^ 

Hence 

m 

Z)(4?+2(*>^ a 3 )) 2 exp(- 7 t) 



fc=i 



max y^(e2 0+1 (t,X, a 2 )) 2 exp(-7t), 
(t r y)e[o,r|xH,J]»- 1 ^ 



fc=l 

Combining this inequality, (I3.20p and the fact that 



^(e 2 fc ?+1 (t,X,6 2 )) 2 exp(- 7 t) < E q+l , 
fe=i 

we deduce 

X^W*'*' «3)) 2 exp(- 7 t) 

k=l 

< exp(-(3(L 2 - S 2 )) max{exp(- p' x S^)E q , exp((3L 2 )E q+1 }, 
Define {3' 2 to be j3[j^ and choose f3 to be — (3' 2 , we infer that 

-/3(L 2 - S 2 ) - ftSi = pS 2 . 

This implies 

m 

^(e^ +2 (t,X,a 3 )) 2 exp(- 7 t) < max{£ fc , E k+l } exp(-/3 2 £ 2 ). (3.21) 
fe=i 

Using the same techniques as the ones that we use to achive (13.201) and 
(13.2 ip . we can prove that 

m 

S^a+j-ifr X ' a i)) 2ex P( - 7^) < max{£ fe , . . . , E k+j _ 2 } exp(- ( 5j_ 1 5 7 _i), 
fc=i 

(3.22) 
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where 

P> = (3[^...^, jG{2,. ..,/-!}. 

Step 2.3: Convergence result 
Setting 

PySi... Sj_i 



2 L 2 . . . 

and 

E k = max {-Efc+j}, 
je{o,...,/-i} L 

we infer from (13351) and ( 13T22D that 

£ fc+ i < £ fc exp(-e),VA; G N. 

Therefore 

-En < -E exp(-ne), Vn G N, 
Hence £4 tends to as k tends to infinity. Which gives that 

m 

lim max V] ||eJ g || c{[ o,T]xR«-ix[a p ,6 p ]) = 0, 

q^oc p=1 j 

and this concludes the proof. 

3.4 Proof of Theorem [3731 

We start proving that 



lim lim / - X 4 | 2 )dt = 0. 



Subtracting (I2.6P and (13.51) . we get 



X t - X, 



t 

q,l 



[ [b(s, X s , 6(s, X,)) - b(s, X!' 1 , 6^(3, X^ l ))]ds (3. 
Jo 

+ f icr((s,X s ,6(s,X s )) - a(s,Xf,6i> l (s,Xf))}dW s , 



o 
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which leads to 



\X t -Xf\ 2 < 2(j\b(s,X s ,9(s,X s ))-b(s,Xf,9«>\ S ,Xf))}ds^ (3.24) 

+2 (fas, X s , 9(s, X.)) - a(s, Xf, 9«> l {s, Xf))}dW?j . 
A simple calculation gives 

E(\X t -Xf\ 2 ) < N 5 E f \b(s,X s ,9(s,X s ))-b(s,Xf,9"' l (s,Xf))\ 2 ds (3.25) 

Jo 

+N 5 E [ \a(s,X s ,9(s,X s )) - a(s,XI' l ,9 q > l (s,Xf))\ 2 ds 
Jo 

< NqE [ [\X s -Xf\ 2 + \9{s,X s )-9«> l {s,Xf)\ 2 ]ds 

Jo 

< NqE [ \X a - Xf\ 2 ds 

Jo 

+ N 7 E [\\9(s,Xf)-9'>> l (s,Xf)\ 2 + \9(s,X s )-9(s,Xf)\ 2 ]ds 
Jo 

< N 8 E [ [\X s -Xf\ 2 + \9{s,Xf)-9«' l {s,Xf)\ 2 ]ds, 

Jo 

where N 5 , iV 6 , N 7 , N$ are positive constants. 

Since {9 q ' 1 } converges uniformly to 9 l , and {9 1 } converges uniformly to 9, 
then for all positive number e, there exists Q{e) such that 

\9(s,Xf)-9«> l (s,Xf)\ < V~e,Vq,l> Q{e). 
Inequality (13.251) leads to 

E(\X t -Xf\ 2 ) < N 8 E [\\X s -Xf\ 2 + e}ds, Vg, I > Q{e). (3.26) 

Jo 

Now, setting 

H{t)= [ E(\X s -X«/\ 2 )ds, 
Jo 

we obtain 

H'(t) < N 9 H(t) + N 9 e, (3.27) 
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where N 9 is a positive constant. Therefore 

H'(t) exp(-N 9 t) - N g H(t) exp(-N 9 t) - N 9 e exp{-N 9 t) < 0. 
This implies 

(H(t) exp(-N 9 t) + eexp(-N 9 t))' < 0, 
and this inequality then leads to 

H(t) exp(-iVgt) + eexp(-A^) < e. 

Consequently, 

H(t) < e(exp(iV 9 T) - 1) 
for q, I greater than Q(e), which leads to 

lim lim / E(\X t - X q ' l \ 2 )dt = 0. 

l^oo q — >oo Jq 

We next prove that 

fT 

lim lim / E(\Y q ' 1 - Y t \ 2 )dt = 0. 



The fact that 



fT 

lim lim / E(\X t - Xf\ 2 )dt = 



implies that Xf' 1 converges to X t almost everywhere. From which, we can in- 
fer that 9 q,l (t, X^' 1 ) converges to 9(t, X t ) almost everywhere. Since {9 q,l (t, Xf' 1 )} 
is bounded, the Lebesgue Dominated Convergence Theorem gives that 



lim lim E [ \Y q > 1 - Y s \ 2 ds = 0, 



and then 



lim lim / E(\Y q ' 1 - Y t \ 2 )dt = 0. 
Now, we prove that 

fT 

lim lim E / \Z q ' 1 - ZA°dt = 0. 



I— >oo q— >oo 
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We get the following equation from the definition of Y t and Y^' , 

Yt-Yf = j\b(s,X s ,Y s )-b(s,Xf,Yf)]ds (3.28) 
+ J\a(s, X s , Y s , Z s ) - a(s, Xf, Y* 1 , Zf)}dW s . 



This implies that 



E J [a(s,X s ,Y s ,Z s )-a(s,X s ,Y s ,Zf)} 2 ds (3.29) 

< AT 10 (E(\Y t - If?) +E j\b(s,X s ,Y s ) - b(s,Xf,Yf)] 2 ds^J 

< Nn (E{\Y t - Y t q ' l \ 2 ) +eJ\\X s - Xf\ 2 + \Y S - Yf\ 2 )ds^j , 



where N w , N u are positive constants. 

Applying Condition (^43) to the inequality (I3.29p . we get 



T 



E J \Z S -Zf\~ds (3.30) 
< N 12 (E(\Y t - Yf\ 2 ) + E J\\X S - Xf\ 2 + \Y S - Y^\ 2 }ds^j , 



where N12 is a positive constant. 
Due to the fact that 



and 



T 

lim lim / E(\X t - Xf\ 2 )dt = 0, 



o 



lim lim E f \Y t - Y t q ' l \ 2 dt = 0, 

l->oo9->oo Jq 



we can see that the RHS of (13.30!) tends to as q and / tend to infinity. Hence 
E /„ \Z t - Zf\«dt tends to as q and / tend to infinity. This concludes the 
proof. 
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4 Conclusion 



We have introduced a new Domain Decomposition Method for a system of 
SDEs. The method has been studied theoretically and proved to be well- 
posed and stable. We have also proposed a new technique to prove the 
convergence of Domain Decomposition Methods for systems of nonlinear 
parabolic equations in n-dimension. The method has the potential to be 
used to prove the convergence of Domain Decomposition Methods for many 
kinds of nonlinear problems. 
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